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o ' 

, Abstract. The results from hydrodynamical TREESPH simulations of 

galaxy clusters are used to investigate the dependence of the final clus- 
, ter X-ray properties upon the numerical resolution and the assumed star 

' formation models for the cooled gas. When cold gas particles are allowed 

^-t , to convert into stars the final gas profiles show a well defined core radius 

If) ' and the temperature profiles are nearly flat. A comparison between runs 

with different star formation methods shows that the results of simula- 
tions, based on star formation methods in which gas conversion into stars 
is controlled by an efficiency parameter c-^, are sensitive to the simula- 
I tion numerical resolution. In this respect star formation methods based 

. instead on a local density threshold, as in Navarro & White (1993), are 

shown to give more stable results. Final X-ray luminosities are found to 

be numerically stable, with uncertainties of a factor ~ 2. 
1 — I I 

o 

:^ 

o 

c/3 . Galaxy clusters are the largest virialized structures known in the universe and 

^ ' also bright X-ray sources. Useful cosmological information can be obtained from 

the statistical properties of the ensemble of X-ray clusters. X-ray observations of 
cluster number counts, the X-ray temperature function (Henry &: Arnoud 1991; 
Edge et al. 1990; Henry 1997) and the X-ray luminosity function (Rosati et al. 
\ 1998; Ebeling et al. 1998) are powerful probes to constrain the values of the 

cosmological parameters ^Iq and erg (Henry & Arnoud 1991; White, Efstathiou, 
& Frenk 1993; Bahcah & Fan 1998; Eke, Cole, & Frenk 1996; Kitayama, Sasaki, 
& Suto 1998) 

Hydrodynamical simulations have been widely used to predict for different 
theoretical frameworks the time evolution of the gas and temperature distribu- 
tions. The numerical methods used are either Eulerian (Cen & Ostriker 1994; 
Anninos & Norman 1996; Bryan &; Norman 1998; Kang et al. 1994; Bryan et al. 
1994, Cen 1997; Cen et al. 1995) with a fixed or adaptative grid, or Lagrangian 
(Evrard 1988; 1990; Thomas & Couchman 1992; Navarro, Frenk, & White 1995; 
Eke, Navarro, & Frenk 1998; Katz & White 1993; Yoshikawa, Itoh, & Suto 
1998; Valdarnini, Ghizzardi, & Bonometto 1999) In these simulations the gas 
component is treated as a single adiabatic fluid, without taking into account the 
effects of radiative cooling, and the physical processes which can be modelled 
are merging, substructure formation, shocks and compressional heating of the 
gas. A comparison between different numerical simulations shows that they are 
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successful in reproducing the gross features of the cluster properties (Prenk et 
al. 1999). With increasing availability of computational power, numerical hy- 
drodynamical simulations have been attempted to model the effects of radiative 
cooling on the gas in the formation and evolution of cluster galaxies (Katz & 
White 1993; Suginohara & Ostriker 1998; Anninos & Norman 1996; Yoshikawa, 
Jing, & Suto 2000; Pearce et al. 2000; Lewis et al. 2000) The numerical prob- 
lems posed by the inclusion of gas cooling are challenging, mainly because the 
required increased spatial resolution also requires that one keeps under control 
two-body heating mechanisms. Previous simulations have produced some con- 
flicting results (Yoshikawa, Jing, & Suto 2000; Pearce et al. 2000), thus the 
question of the minimum resolution in this kind of simulations is still to be fully 
settled. 

In this contribution I will show the preliminary results that have been ob- 
tained from a series of hydro dynamical SPH simulations of galaxy clusters. The 
simulations have different numerical resolution and include the effects on the 
gas component of radiative losses, star formation and energy feedback from SN. 
Final profiles arc compared in order to assess the effects of numerical resolution, 
or different star formation prescriptions, on the cluster X-ray variables. 

2. Simulations 

In a previous paper (Valdarnini et al. 1999, hereafter VGB) a large set of 
hydrodynamical simulations was used to study the global X-ray cluster mor- 
phology and its evolution. The simulations were run using a TREESPH code 
with no gas cooling or heating. In order to assess the numerical reliability of 
the numerical integrations, four different clusters were selected as a represen- 
tative sample of all the simulation clusters. For this cluster sample a large set 
of different integrations was performed by varying two numerical input param- 
eters: the number of particles and the softening parameters. I refer to VGB 
for a detailed description of the simulations. In Valdarnini (2001) I have used 
the same cluster sample to study the effects of including in the simulations ad- 
ditional physics such as gas cooling and star formation. I will report here the 
results for the cluster with label ACDMOO in VGB. This is the most massive 
cluster (M„ ~ 1.5- lO^'^^AfQ) extracted from a cosmological ACDM N-body simu- 
lation with size L = 200h~^ Mpc, matter density = 0.3 and Hubble constant 
Hq = 7QKm sec^^ Mpc~^. In order to check the effects of radiative cooling for 
this cluster a set of TREESPH simulations was performed, with initial condi- 
tions provided by the cosmological simulation, and different values of the gas 
softening parameter Eg and number of gas particles Ng. For this cluster Table 1 
reports the values of £g and Ng for the simulation runs. The generic simulation 
has cluster index clOO — j, with j = 00,01, ...,05. The cluster clOO — 00 is the 
reference case without cooling. The effects of radiative cooling are modelled in 
these simulations by adding to the SPH thermal energy equation an energy-sink 
term. The total cooling function includes contribution from recombination and 
collisional excitation, bremsstrahlung and inverse Compton cooling. 

Allowing the gas to cool radiatively will produce dense clumps of gas at low 
temperatures (~ 10^ °K). In these regions the gas will be thermally unstable and 
will likely meet the physical conditions to form stars. In TREESPH simulations 
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Figure 1. Final density and temperature profiles in the simula- 
tion runs including radiative cooling and star formation. Density is in 
units of the critical density. Different curves are for integrations with 
different numbers of particles and different softening parameters (see 
Tabic 1). The simulation run with index —00 is the integration without 
cooling. 



star formation (SF) processes have been implemented using different approaches. 
According to Katz, Weinberg, & Hernquist (1996, hereafter KWH) a gas particle 
is in a star forming region if the flow is convergent and the local sound crossing 
time is larger than the dynamical time (i.e. the fluid is Jeans unstable). In a 
simplified version Navarro & White (1993 , hereafter NW) assume as a sufficient 
condition that a gas particle must be in a convergent flow and its density exceeds 
a threshold {pg > pg^c = 7 • lO^^^^rcm"^). 

If a gas particle meets these criteria then it is selected as an eligible particle 
to form stars. The local star formation rate (SFR) obeys the equation 

dpg/dt = -C^pg/Tg = -dp^jdt , (l) 

where pg is the gas density, is the star density , is a characteristic dimen- 
sionless efficiency parameter, Tg is the local collapse time and it is the maximum 
of the local cooling time Tc and the dynamical time r^. Gas particles with 
T ^ W^°K have long cooling times and Tg = t^. At each step the probability 
that a gas particle will form stars in a time step Ai is compared with a uniform 
random number. If the test is successfull then a mass fraction of the gas mass 
is converted into a new collisionlcss particle. This star particle has the position, 
velocity and gravitational softening of the original gas particle. Typical assumed 
values are = 1/3 and = 0.1 (KWH). For the local SFR, NW adopt Eq. 1 
with Ci, = 1, Tg = Td and 8^, = 1/2 when a gas particle can convert part of its 
mass into a star particle. 
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The numerical tests have been performed following the NW prescriptions 
for selecting gas particles which can form stars. Once a star particle is created it 
can release energy into the surrounding gas through supernova (SN) explosions. 
All the stars with mass above 8Mq end as a SN, leaving a IAMq remnant. 
The SN energy (~ 10^^ erg) is released gradually into the gas according to the 
lifetime of stars of different masses. 



3. Results 

The final radial gas density and temperature profiles arc shown in Fig. 1. As 
can be inferred from Table 1, the numerical strategy has been to perform runs of 
increasing spatial resolution in order to resolve the core radius of the gas density 
profile. The most important result is that the inclusion of a star formation 
prescription has been effective in removing the cold gas particles (T ^ 10^ °K) 
from the cluster center. In all the simulations the gas density profiles have a 
well-defined core radius, with size rc — 50 — lOOKpc, approximately 0.05 of 
the virial radius (~ 2.9Mpc). The shapes of the temperature profiles show 
that convergence is achieved for Ng ^ 20, 000. All the central values for the gas 
temperatures at r = lOKpc are within a factor ~ 1.5 for the highest resolution 
simulation runs. The profiles increase inwards from the virial radius up to ~ 
100 — 200Kpc. Thereafter the profiles stay almost flat, or with a modest decline 
in T(r) towards r = 0. There is not a strong drop of the temperature in the very 
central region. Final X-ray luminosities are quite stable {Lx — 2- 10"^^ ergsec~^) 
versus the simulation numerical resolution. Compared to the non-radiative run 
the luminosities increase by a factor ~ 2. 

A different question is the sensitivity of the estimated cluster properties to 
the numerical method used to describe star formation in the hydrodynamical 
simulations. To this end three simulations were performed with the same nu- 
merical parameters as for the clOO — 05 run, but using different SF methods or 
parameters. The NW simulation (I) is the standard case showed in Fig. 1. In 
the other two runs (II and III) the KWH prescription is adopted for converting 
gas particles into stars, but with different values of the star formation efficiency 
parameter c^: 0.1 and 1. (see Fig. 2). The most important differences are found 
between the KWH simulations with different c^,. The differences are dramatic 
in the final X-ray luminosities, which differ by a factor ~ 40. The source of this 
discrepancy relies in the different gas density profiles, which have substantial 
differences in the cluster core regions for r ^ lOOKpc. These differences are lo- 
calized at the cluster center, beyond r ~ lOOKpc all the profiles converge, as it 
is shown in the plots of Fig. 2. The temperature profiles have a peak value of 
2± 10^ °fC at r ~ lOOKpc and thereafter decline outward by a factor ~ 2 out to 
ry. Below ^ lOOKpc the profiles instead show large differences. Compared to 
the NW run the simulation with = 0.1 has gas temperatures which decline 
inwards by a factor ~ 10 from ~ lOOKpc down to r ~ lOKpc. These radial 
decays follow because of the less efficient conversion of the cooled gas into stars 
compared to the NW run. There is a remarkable agreement for the ratio of 
cluster mass locked into stars to the gas mass, which is ~ 10% at for all 
the runs considered. In order to assess the effects of numerical resolution upon 
final results simulations I, II and III have been run again but with a number 
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Figure 2. Plots showing several cluster properties in simulation runs 
with different SF prescriptions. The simulation parameters are those of 
clOO — 05 (see Table 1). The continuous line refers to the NW method, 
the others to KWH with different ( c in the bottom left panel) . Top 
left: ratio of the star mass within the radius r over the gas mass within 
r. Top right: final radial density behavior for the gas component. The 
simulation results are compared with the corresponding high-resolution 
runs (H). The H simulations have the same SF parameters of the par- 
ent simulations, but the numerical parameters are given in the last row 
of Table 1. To facilitate a comparison the radial profiles of the high 
resolution results have been shifted downward by lO'^, k = 2 for den- 
sities and k = 1 for temperatures. Bottom right: Radial temperature 
profiles. Bottom left : X-ray luminosity versus time, the thick lines 
correspond to the high resolution simulations. 
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of particles increased by a factor ~ 3 (clOO — 05H). These simulations will be 
referred as IH, IIH and IIIH, respectively. The simulation results are shown in 
Fig. 2. For simulations IH there are not appreciable differences in the radial 
profiles. This confirms the previous results, that is the NW runs have reached 
numerical convergence in the physical variables for the numerical parameters of 
the clOO — 05 simulation of Table 1. The profiles of simulation IIH are instead 
different from those of run II at r ;S 50Kpc. The strong drop in T(r) has been 
removed and the gas density profile is much closer to the NW one. Simulation 
IIIH yields final profiles very similar to the ones of the parent simulation. The 
bottom left panel of Fig. 2 shows that high resolution runs have final X-ray 
luminosities which can differ within a factor ~ 2 from the parent simulations. 

To summarize, the above results demonstrate that simulations I and III have 
an adequate numerical resolution to reliably predict X-ray cluster properties, 
such as the X-ray luminosity. For simulation II ( KWH with = 0.1 ) there are 
large differences at the cluster core between the final profiles when the numerical 
resolution is increased. For the simulation runs I and HI Lx — 10^^ ergsec~^, 
while Lx — A-W^^ergsec^^ for the run with c^, = 0.1. Thus simulations I and III 
are consistent at the la level with the — Tx relation derived from a sample 
of cooling flow clusters (Allen & Fabian 1998). For simulation II the predicted 
cluster temperature is outside the 2a limits. 
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Table 1. Simulation parameters of the test runs for the ACDMOO 
cluster. clOO — 00 is the reference case with no cooling, taken from 
VGB. gravitational softening parameter for the gas in Kpc. ^: 
mass of the gas particles in h~^MQ (the cosmology is for = 0.3 and 
h = 0.7). : mass of the dark particles. number of gas particles. 
^ : as in the previous column but for dark particles. total number 
of simulation particles. value of the treecode gravitational tolerance 
parameter. ^ : initial redshift for the simulation. The last row gives the 
numerical parameters for the high-resolution runs used to test different 
SF methods. For these runs gravitational quadrupole corrections were 
enabled. 
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